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ABSTRACT 

This  memorandum  describes  a  model  and  computer  program  designed  to 
compute  kill  probabilities  for  certain  firing  tactics  waged  against  an  evading 
submarine.  The  model  was  developed  assuming  that  the  attack  is  imminent  and 
that  the  weapon  will  be  directed  at  the  point  where  the  submarine  was  last  con¬ 
tacted.  The  model's  design  is  centered  around  a  determination  of  sets  of  points 
which  represent  the  locus  of  the  evading  submarine.  The  locus  is  determined  by 
variation  of  a  simple  evasion  tactic.  Theaimpoiiit  is  considered  to  be  circularly, 
normally  distributed  about  the  true  position  of  the  submarine  at  the  time  of  last 
contact. 
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1.0  Discussion  of  the  Model: 

This  model  and  program  were  developed  to  provide  a  means  for  computing 
kill  probabilities  for  certain  ASW  tactics  .  The  model  and  program,  though 
predicated  upon  a  particular  tactical  situation,  are  sufficiently  flexible  to  handle 
a  variety  of  cases  and  should  permit  design  modifications  without  major  revision. 

1. 1  Tactical  Situation: 

The  tactical  situation  which  this  model  is  designed  to  treat  is  described 
as  follows: 

A  helicopter,  assisting  a  delivery  aircraft  in  an  attack,  has  a  firmly 
established  sonar  contact  with  a  submarine. 

1.1.0  Condition  1: 

If  it  is  found  that  the  assisting  helicopter's  distance  from  the  sub  is  within 
a  minimum  safe  distance,  certain  tasks  must  be  performed  before  the  assisting 
helicopter  can  retire,  under  emergency  power,  to  the  minimum  safe  distance. 

The  delivery  aircraft  must  await  communications  and  directions  from  the 
assisting  helicopter,  and  then  delay  for  at  least  some  maximum  time  before 
maneuvering  to  the  predicted  position  and  dropping  the  weapon. 

1.1.1  Condition  2: 

If  it  is  found  that  the  assisting  helicopter  is  outside  the  minimum  safe 
distance  from  the  submarine's  position,  the  assisting  helicopter  tracks  the 
submarine  and  continuously  vectors  the  delivery  aircraft  until  the  weapon 
is  dropped. 

1.2  Analysis: 

The  key  factor  in  these  situations  is  the  time  that  the  submarine  remains 
out  of  the  sonar  contact.  That  time,  under  condition  2,  is,  usually,  negligibly 
small,  even  if  a  delay  maneuver  is  incorporated.  Not  many  escapes  are 
possible  here.  The  situation  defined  by  condition  1  is  such  that  a  submarine's 
chance  to  escape  is  measurably  enhanced.  The  probability  of  escape  is  a 
function  of  the  weapon's  lethal  radius  and  the  distance  between  the  submarine 
and  the  weapon's  burst  point  at  detonation. 

The  locations  attained  by  a  submarine,  for  all  evasion  tactics  executed 
within  a  given  time  interval,  will  be  found  at  various  points,  more  frequently 
in  one  spot,  less  frequently  in  another,  at  a  greater  or  lesser  distance  from 
the  origin,  taken  as  the  point  of  last  sonar  contact.  The  problem,  then,  is 
to  find  the  values  for  the  kill  probability  for  distributions  of  points  in  the  locus 
of  the  evading  submarine. 


2.  0  Method  of  Solution: 

The  program  computes  sets  of  position  which  lie  in  a  lobe  of  the  locus  of 
points  of  the  evading  submarine.  The  origin  of  the  locus  is  the  point  of  last 
contact,  and  the  sub  is  assumed  to  be  initially  moving  along  the  positive  y-axis. 
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The  terminal  positions  of  the  submarine  are  found  from  the  initial  speed  of  the 
submarine  at  the  beginning  of  the  blind  time,  the  blind  time  interval,  andvar- 
iations  of  a  simple  evasion  tactic.  The  blind  time,  TB,  is  the  time  between  the 
severance  of  sonar  contact  and  the  weapon's  detonation.  The  terminal  positions 
so  found  comprise  a  representative  sample  of  the  distribution  of  points  in  the 
locus  of  an  evading  submarine.  The  evasion  tactics  are  described  as  follows: 

(a)  During  a  multiple,  t,  of  a  fraction  of  the  blind  time,  TB,  the  submarine 
cruises  on  course  at  speed,  V0  .  The  value,  (t),  is  reckoned  from  the  beginning 
of  the  blind  time,  TB. 

(b)  After  time,  t,  the  submarine  executes  a  turn  with  radius  of  turn,  R, 
through  a  multiple,  A,  of  a  fraction  of  the  turn  angle,  0  ,  at  an  average  speed, 

V,  a  function  of  V0  and  A  . 

(c)  If  there  is  any  remaining  blind  time,  the  submarine  cruises  along  its 
new  course  at  speed,  V0,  until  all  of  the  blind  time  is  consumed. 

(d)  Then,  for  the  given  combination  of  TB,  V0,  t  and  A,  a  terminal  po¬ 
sition  is  established  for  the  submarine. 

(e)  The  process  is  repeated  for  all  combinations  of  a  finite  set  of  values 
in  the  range  of  t  and  a  finite  set  of  values  in  the  range  of  A,  where  0<  t  sTB 
is  the  range  of  t  and  0£  As  0  is  the  range  of  A  . 

2. 1  The  equations  for  calculating  the  coordinates,  (X^.  ^  ,  Y^.  ^  ),  of  the  po¬ 
sitions  in  the  locus  of  the  evading  target  are  derived  as  follows:  (see  figure  1, 
appendix  D) 

Given: 

TB  =  blind  time 
e  =  turn  angle 
M  =  a  denominator 

AT=  TB/M,  a  fraction  of  the  blind  time 

N.  =  an  integer  which  is  allowed  to  take  on  successive  values  in 
the  range,  0  s  N.  £  M 

t  =  N.  AT,  a  multiple  of  a  fraction  of  the  blind  time 

L  =  a  denominator 

A6  =  0  /L,  a  fraction  of  the  turn  angle 

P.  =  an  integer  which  is  allowed  to  take  on  successive  values  in  the 
J  range,  0  s  P .  s  L 

A  =  Pj  A  0 ,  a  multiple  of  a  fraction  of  the  turn  angle 
R  -  turn  radius 
V0  =  initial  speed  of  submarine 
.  V  =  average  speed  of  submarine  during  turn 
Vj=  the  instantaneous  speed  of  a  submarine  during  a  turn 
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First,  let  us  derive  an  expression  for  Vj.  Assume  that  the  instantaneous  speed 

of  a  submarine  during  a  turn  is  a  continuous,  linear  function  of  the  turn  angle,  9, 
and  that  the  speed  decreases  to  55  percent  of  its  initial  value  at  0  =  90°.  These  as¬ 
sumptions  are  consistent  with  the  results  of  studies  made  on  submarine  turning 
rates.  Then, 


(V0,  0)  and  (0. 55V0,  tr./2)  are  points  on  the  line  and 
(Vj-VoVe  =  —2 V0  (1—0.  55)/tt  =  — 0.90Vo/tt 

ttVj  —  TiVo  =  — 0.90Vo  0 

Vj  =  (nVo  -  0.90Vo  0)tt  =  -V0(l  -  0.90  0/n) 

Then,  the  average  value,  V,  of  the  function,  VT,  between  0  and  0  is, 
0:  1 

V  =  £  Vo  (1-0.900/rr)  d0 

V=  Vo/0  |  0-O.9O02/2tt  I  o  =(Vo/0)  .  (0  -o.902/2n) 


Now, 


Then, 


V  =  Vo  (1-O.9O0/2tt)  . 


TO  =  R(P^.A9)/V  =  the  time  to  turn  through  the  angle,  A  =  PA8 
TI  =  TB-NjAT  =  the  time  remaining  after  the  interval,  t  =  RAT 


T2  =  Tl-TO  =  time  remaining  to  cruise  on  new  course. 


X(i,  j)  ~  R  (1  “  cos  (VT1/R)) 
Y(i  =  Vp N.A T  +  R  sin  ( VT1/R) 

X(.  R  (1— cos  (P^A0)) 

Y(i  j>  =  VoN.AT  +  Rsind’jA0) 


JT2  <0 


^T2  =0 


(1) 


(2) 
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X  ^  =  R  (1-cos  (Pj  A  0 ))  +  V0T  2  sin  (P  A  0) 

Y  -  V0N.  A T  +  R  sin (P.  A  0)  +  V0T  2  cos  (P.  A  0) 


»  T2  >  0 
y 


(3) 


Equation  (1)  is  used  if  the  turn  cannot  be  completed;  equation  (2)  if  it  can  but, 
in  so  doing,  consumes  all  of  the  blind  time;  and  equation  (3)  if  all  aspects  can  be 
accomplished  within  the  blind  time. 


2.  2  There  now  remains  the  problem  of  determining  the  point  at  which  the  attack 
on  the  target  is  aimed.  That  point  is  the  point  at  which  the  target  was  last  observed. 
The  observations  of  that  point  are  assumed  to  be  circularly,  normally  distributed 
about  the  true  position  of  the  submarine  at  T  =  0,  the  beginning  of  the  blind  time. 
Coordinates  of  the  aimpoint  can  be  established  as  follows: 


Assume  a  polar,,  coordinate  system  with  origin  at  the  sighting  ship,  at  T  =  0, 
and  the  polar  coordinates  of  the  target  are  (R^,  tt/2).  (See  figure  2),  appendix  D.  Then, 


(a) 

<b) 

(c) 

(d) 

(e) 

(f) 

(g) 


(Rt  ,  tt/2)  =  polar  coordinates  of  the  target's  true  position. 

(Ro  .  oto)  -  polar  coordinates  of  the  target's  observed  position 
AR  =  Ro  -  Rt  =  error  in  the  range  observation 

Aa  =  ao  -  tt/2  =  error  in  the  azimuth  observation  then,  transforming 
to  cartesian  coordinates, 

X£  =  Rt  cos  tt/2 

Yt  =  Rt  sin  tt/2  =  R^ 

Xo  =  Ro  cos  a0 
Yo  =  Ro  sin  oo 


}  cartesian  coordinates  of  the 
target's  true  position 


Vcartesian  coordinates  of  the 
(target's  observed  position 


but,  from  item  (d), 

(h)  ao  =  tt/2  +  Aa,  and  therefore 

(i)  Xo  =  R0  cos  (rr/2  +  Aa)  =  R0  sin  Aa 

(j)  Y0  =  R0  sin  (tt/2  +  Aa)  =  R0  cos  Aa 


then,  by  translation  of  axes, 

(k)  Xo  =  Xo  -  X£  =  R0  sin  Aa 

(l)  Yi  =  Y0  -  Yt  =  Ro  cos  Aa  -  Rt 

then,  (Xo  ,  Yp  )  =  coordinates  of  the  aimpoint . 
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The  coordinates  of  the  aimpoint  with  origin  at  the  sub  at  T  =  0,  are  found  for 
each  variation  of  the  evasion;  they  are  computed  as  follows: 

Given, 

cr  =  standard  deviation  in  range  error  (2%  of  range) 

K 

cr  =  standard  deviation  in  azimuth  error 
a 

then,  choose 

G^  =  a  random  number  from  a  Gaussian  distribution  of  random  numbers 
with  mean  zero  and  standard  deviation  equal  to  1. 

and  choose 

G  =  a  random  number  from  a  Gaussian  distribution  of  random  numbers 
a  with  mean  zero  and  standard  deviation  equal  to  1. 


then, 


Rt  =  crR  .  Gr/  0.02  =  50  .  a-R  .  G^  -  true  range  of  aimpoint  from 
origin  at  assisting  ship  at  T  =  0. 


Gn  =  an  observed  range  from  origin  at  assisting 


ship  at  T  =  0. 


and, 


A  a  =  cr  .  G  =  the  azimuth  error 
a  a 


(Xi,  Y o)  =  (-Ro  sin  Aa,  R0  cos  Aa  -  R  )■. 

Then,  the  distance  between  the  target's  position  at  T  =  TB  and  the  aim 
point  is: 


-x&2  + 


(Y 


(i,  j) 


2.3  This  program  uses  the  SHARE  subroutine  OCIP.  OCIP  uses  the  criterion 
of  circular  coverage  to  establish  the  probability  that  a  Weapon  with  a  certain 
radius  of  effectiveness  and  a  ballistic  dispersion  with  standard  deviation  of  crD 

D 

will  effectively  dispose  of  a  target  randomly  displaced  from  the  aimpoint.  The 
distribution  of,  the  target's  displacement  is  characterized  by  a  circular,  normal 
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distribution  of  standard  deviation  o-p  .A  The  ballistic  dispersion  and  the  target's 
displacement  are  centered  at  the  aimpoint.  For  our  purposes,  D^.  .j  is  the  target's 
displacement  and  the  weapon's  lethal  radius,  RKILL,  is  the  radius  of  effectiveness. 

Input  to  OCIP  must  be  given  in  the  form: 


(D(.  j)/(r,  RKILL/o-,  PK), 


where, 


-l 


2  .  2 
+<TP 


-  <ru 


3 


■*) 


and,  the  output  parameter,  PK,  in  our  usage,  is  the  probability  that  the  firing 
tactic  results  in  a  kill. 

2.  4  Consequently,  D^.  ^/<r  and  RKILL/c  can  be  found  for  each  position  in  the 

sample  of  computed  positions  in  the  locus  of  The  evading  target  and,  using  the 
circular  coverage  function  as  programmed  in  OCIP,  a  kill  probability  for  the 
individuals  in  the  sample  can  be  calculated.  The  average  of  the  kill  probabilities 
of  the  individuals  in  the  sample  is  taken  as  the  kill  probability  for  the  tactic 
against  an  evading  target.  Kill  probabilities  are  computed  for  several  values  of 
lethal  radius  and  submarine  speed. 


2.  5  This  model  assumes  that  the  time  to  turn  is  uniformly  distributed  throughout 
the  blind  time  and  that  the  turn  angle  is,  also,  uniformly  distributed  through  all 
possible  values  of  that  angle.  The  model  was  constructed  to  calculate  the  kill 
probability  regardless  of  the  evasion  tactic.  Consequently,  the  evasion  tactic 
contained  in  the  model  is  designed  only  as  a  device  for  generating  sets  of  points 
which  define  a  set  of  limiting  distances  attained  by  a  target  in  a  given  time, 
without  giving  special  regard  to  favored  evasions.  The  blind  time  is  expected 
to  be  of  short  duration.  The  blind  time,  thus,  does  limit  the  design  of  a  more 
sophisticated  evasion  which  one  could  expect  a  submarine  to  accomplish  under 
imminent  attack. 

A  study  of  the  results,  oriented  toward  studying  evasion  tactics,  should 
help  to  improve  the  model  and  may  well  suggest  better  countermeasures  against 
evasive  maneuvers. 

Consideration  is  currently  being  given  to  design  modifications  which  will 
enable  specification  of  favored  evasion  schemes  and  provide  for  prediction  of 
aimpoints  involving  bearing  and  speed  prediction  errors,  offset  in  range  and  azimuth. 

3.0  Program  Input: 

The  program  will  accept  card  input.  At  least  two  cards  are  required,  a  control 
card  and  a  test  case  card.  A  control  card  can  be  followed  by  any  number  of  test 
cases.  Any  number  of  test  sets  may  be  submitted. 
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3. 1  The  Control  Card: 

A  set  of  test  case  cards  must  be  preceded  by  a  control  card,  The  control  card 
is  used  to  indicate  the  aumbe*r  of  sets  and  the  number  of  test  eases  per  set  in  the 
input  deck,  and  the  output  options.  There  are  three  output  options.  See  paragraph  5.0, 
Program  Output. 

3 . 1. 0  The  Test  Case  Card; 

There  are  eleven  input  parameters.  The  eleven  parameters  are  described 
as  follows: 

3.1.1  Maximum  Turn  Angle; 

This  is  the  maximum  angle,  Q°,  through  which  the  evasive  turn  will  be  executed, 
A  fraction  of  this  angle  is  taken  and  successive  multiples  of  that  fraction  are  used 
to  interate  the  computation  of  the  kill  probability.  Theta  must  be  given  in  degrees. 

3.1.2  Multiples  of  Turn  Angle: 

This  parameter,  L,  is  used  by  the  program  to  determine  the  number  of  values 
of  Pj  Aq  through  which  the  submar  ine  will  be  programmed  to  turn.  A0  =  0°/L  and 

0  s  Pj  s  L. 

3.1.3  Blind  Time: 

Blind  time,  TB,  is  the  length  of  time  between  the  bomb  burst  and  the  time  when 
sonar  contact  is  severed.  A  fraction  of  this  value  is  determined,  using  M  (see 
paragraph  3.1.4),  and  successive  multiples  of  that  fraction  are  used  to  iterate  the 
computation  of  the  kill  probability.  Blind  time  is  given  in  seconds, 

3.1.4  Multiples  of  Time: 

This  parameter,  M,  is  used  by  the  program  to  determine  the  number  of  values 
of  N.  AT  so  that  the  moment  in  time  at  which  the  evasive  turn  will  be  executed  will 

be  uniformly  distributed.  AT  =  TB/M  and  0  s  N.  s  M. 

3.1.5  The  Turn  Radius: 

The  turn  radius,  R,  is  the  radius  of  the  evasive  turn.  This  value  is  given 
in  yards. 

3.1.6  Sub  Speed: 

This  is  the  speed  of  the  submarine.  The  speed,  V,  is  given  in  knots.  Five 
successive  multiples  of  V/5  are  used  as  the  sub  speed  by  the  program. 

3.1.7  The  Damage  Radius: 

(RKILL)  This  specifies  the  distance  from  the  point  of  the  burst  of  the  bomb 
within  which  a  1.00  percent  probability  of  kill  will  be  satisfied.  RKILL  is  given 
in  feet.  Ten  successive  multiples  of  RKILL/ 10  are  used  as  the  lethal  radius  by 
the  program. 
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3.1.8  The  Sonar  Azimuth  Error: 

<r  is  given  in'  degrees  and  is  the  standard  deviation  of  the  error  in 

3. 

observation  due  to  azimuth. 

3.1.9  The  Sonar  Range  Error: 

cr^  is  given  in  feet.  It  is  the  standard  deviation  of  the  error  in 
observation  due  to  range. 

3.1.10  The  Ballistic  Error: 

orD  is  given  in  feet  and  is  the  standard  deviation  of  the  ballistic  error. 

D 

3.1.11  Displacement  Error: 

o-p  is  given  in  feet  and  is  the  standard  deviation  of  the  circular,  normal 
distribution  of  the  target's  displacement  from  the  aimpoint.. 


4.0  Key  Punching  Instructions: 


4. 1  Control  Card: 

There  are  four  fields  on  a  control  card.  The  first  three  fields  are 
six  columns  long,  and  all  entries  therein  must  be  right-justified.  The 
fourth  field  is  fifty-four  columns  long  and  may  be  used  for  remarks. 

Cols.  1-6:  This  field  is  used  to  specify  the  number  of  test  case 
cards  under  the  purview  of  the  given  control  card. 

Cols.  7-12:  This  field,  if  blank  or  punched  zero,  indicates  that 

no  other  test  sets  follow  the  given  set;  if  punched  with 
a  one,  another  test  set  is  expected  to  follow.,  A  test 
set  consists  of  a  control  card  followed  by  any  number 
of  test  case  cards. 


Cols.  13-18:  This  field  is  used  to  specify  the  output  options;  if  blank 
or  zero,  the  average,  summary  results  of  the  compu¬ 
tations  from  the  N.  x  combinations  of  the  multiples  of 

blind  time  and  turn  angle  will  be  printed;  if  punched  with 
a  one,  the  individual  results  from  the  computation  on  the 
N.  x  Pj  combinations  of  blind  time  and  turn  angle  will  be 

printed  in  addition  to  the  average,  summary  results;  if 
punched  with  a  two,  the  zero-time  combinations  of  blind 
time  and  turn  angles  will  be  printed  in  addition  to  the 
average,  summary  results. 
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4. 2  Test  Case  Card: 

There  are  eleven  fields  on  a  test  case  card.  Each  field  is  six  columns  long, 
and  any  entry  in  a  field  must  be  right- justified.  The  card  should  be  punched  as 
follows: 

Max  Turn  Angle  (0°  ),  degrees. 

Multiples  of  0U  (L). 

Blind  Time  (TB),  seconds. 

Multiples  of  Blind  Time  (M). 

Turn  Radius  (R),  yards. 

Sub  Speed  (V),  knots. 

Damage  Radius  (RKILL),  feet. 

Azimuth  Error  (o^),  degrees. 


Cols. 

1-6: 

Cols. 

7-12: 

Cols. 

13-18: 

Cols. 

19-24: 

Cols . 

25-30: 

Cols. 

31-36: 

Cols, 

37-42: 

Cols. 

43-48: 

Cols, 

49-54: 

Cols. 

55-60: 

Cols. 

61-66: 

Sonar  Range  Error  (<Tg),  feet. 
Ballistic  Error  (cr0),  feet. 

D 

Displacement  Error  (<Tp),  feet. 


5.0  Program  Output: 

There  are  three  output  options.  These  options  are  set  by  key-punching 
the  PBIT  field  (cols.  13-18)  on  the  control  card.  The  output  is  designed  so 
that  kill  probabilities  can  be  read  out  over  ranges  of  speed  and  lethal  radius. 


5.1  Option  A:  (PBIT  field  is  blank  or  zero) 

The  following  are  printed  under  option  A. 

(a)  a  value  of  the  sub's  speed  in  knots  equal  to  k  x  V/5,,  where 
1  i  k  s  5  and  V  is  the  value  of  the  speed  given  in  the  input. 

(b)  a  value  of  lethal  radius  in  feet  equal  to  1  x  RKILL/ 10,  where 
1.0s  1  s  10  and  RKILL  is  the  damage  radius  given  in  the  input. 

Pi  2 

(c)  sigma,  equal  to  y/o-g  +  o-p 

(d)  for  each  of  the  k  x  1  combinations  of  speed  and  RKILL,  an  average 
kill  probability  (PK),  computed  from  the  combinations  of  N.  time 

intervals  and  P^  angle  intervals  is  printed. 

(e)  the  total  of  the  number  of  individual  computations  of  PK  involved  in 
computing  the  average.  Total  =  (M  + 1. 0)  x  (L+ 1.0). 
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5. 2  Option  B:  (PBIT  field  is  set  to  one) 

For  each  of  the  k  x  1  combinations  of  speed  and  RKILL,  the  following  will 
be  printed  under  option  B  in  addition  totheoutput  in  option  A: 

(a)  the  time  intervals;  (M  +  1.0)  values. 

(b)  the  coordinates  of  the  position  of  the  submarine  from  the  origin  in 
feet;  (M  +  1.0)  x  (L  +  1.0)  values. 

(c)  the  distance  of  the  submarine  from  the  origin  in  feet; 

(M  +  1 . 0)  x  (L  +  1 .  0)  values. 

/2  2 

(Tg  +  <Tp  . 

(e)  the  fraction  of  the  input  turn  angle;  (L  +  1.0)  values. 

(f)  the  fraction  of  the  input  speed;  k  values. 

(g)  the  fraction  of  the  input  lethal  radius;  1,  values. 

(h)  the  kill  probabilities,  (M  +  1. 0)  x  (L  +  1.0). 

5.3  Option  C;  (PBIT  field  is  set  to  two) 

Option  C  is  the  same  as  option  B  but  involves  only  those  results  concerning 
the  zero  time  interval;  that  is  item,  a,  under  option  B  is  always  zero  in  option  C 


6.0  User's  Instructions: 

A  user  must  submit  the  following  information  for  each  test  set: 

(a)  the  number  of  cases  in  the  test  set. 

(b)  indication  whether  or  not  the  given  set  is  the  last  set  in  the  input  deck. 

(c)  indications  of  the  output  option  desired. 

(d)  for  each  test  case,  values  for:  turn  angle  in  degrees,  multiples  of 
turn  angle,  blind  time  in  seconds,  multiples  of  blind  time,  turn  radius 
in  yards,  speed  of  sub  in  knots,  lethal  radius  in  feet,  azimuth  error  in 
degrees,  range  error  in  feet,  ballistic  error  in  feet,  and  displacement 
error  in  feet. 

This  information  can  be  conveniently  given  on  an  agendum  sheet  such  as 
that  shown  in  appendix  B.  The  output  shown  in  appendix  C  was  derived  from 
the  sample  input  given  in  appendix  B.  In  addition  to  the  circular  coverage  sub¬ 
routine,  OCIP,  the  program  uses  the  random  number  generator,  GRNUMB. 

The  binary  decks  for  these  subroutines  must  be  incorporated  with  the  binary 
deck  for  the  program  described  in  this  document. 


12 


INTERIM  RESEARCH  MEMORANDUM 

IRM-25 


APPENDIX  A 
FORTRAN  LISTING 
AND 

FLOW  DIAGRAM 
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C  COMPUTE  ASH  KILL  PROBABILTY 

DIMENSION  REMARK! 9) 

1  FORMAT  1216, F6.0»9A6) 

3  FORMAT! 1H  ,  42H  COMPUTATION  OF  EFFECTIVE  KILL  PROB AB I L I T Y// IH  , 

L50X,  19HIINPUT,  PARAMETER  SET/1H  IOHTURN  ANGLE  ,  1 X  .  9HMULT IPLES  ,  2X,  5HBL 

2  IND6X9HMULT  IPLES2X4HTURN7X3HSU88X6HDAMAGE5X7HAZ I MUTH4X5HRANGE6X9HB 

3  ALL  I  ST IC2X6HT  ARGE  T/  IH  ,  8H  THETA  ,5X,8H0F  THETA,  3X,  4HT  IME  ,.5X,  7HOF 

4  T I ME4X, 6HR AO  1  US , 9X , 5HSPEED , 4X , 6HR  AD  I  US , 6X , 5HERROR , 6X , 5HERROR, 8X ,5 
5HERROR , 4X , 5HERRQR ) 

4  FORMAT ( 1 1F6 . 0  ) 

5  FORMAT!  IH  10  I  F 10.  2  ,.IX  ) ,  F 8. 2) 

6  FORMAT  (  IH  ,  3X  ,  4HT  IME  ,.6X  ,.5HTHE  TA  ,.4X  ,  5HXSUB  1 , 5X  ,  5HYSUBI  ,  5X,  5HDSUBI, 
15X,5HSIGMA,5X,5HRKILL,SX,5HSPEED,7X,2HPK) 

7  FORMAT ( 1H0, 30X, 16HC0MPUTAT ION  SETS) 

21  FORMAT! IH1, 216, F6.0.9A6) 

PL-3., 1415927 

22  READ  l  (NCASES,  NSETS,  P8I  T,  I  REMARK I  I  )  ,  1—1,9) 
PRINT2I,NCASES,NSETS,PBIT,'(REMARK(  1  ),  1*1,9) 

2  READ  4, ARC, ARCN,  BT ,  BTM  .  ARCRAD  ,  VSUB  ,  DAMRAD ,  S  A,, SR ,  SB ,  SP 

C  PRINT  HEAD  I NG  AND  I.NPUT  SET 

PRINT  3 

PRINTS, ARC , ARCN, BT , BTM, ARCRAD , VSUB , DAMRAD , SA , SR , SB , SP 
PRINT  7 

LF I PB IT- 1 .0 ) 56, 55 , 55 

55  PRINT  6 

56  C-15.0 
NARG— ARCN *1.0 
MBT-BTM* 1 . 0 
VSU8L=1.689*VSUB/5.0 
SA=PL*SA/180.0 
TRAO- 3.0* ARCRAD 
ALPHA=PI*ARC/  (  180-  0*ARCN ) 

OT-B  T/8TM 

S=.SQKTF(SB**2+SP*»2) 

2-1.0 

00  153  L-1,5 

FVSUB=Z*VSU8I 

SPEE0-FVSUB/1.689 

2=2*1. Q 

R-1.0 

00  152  M-1,10 

8  SUM-0.0 
SUMPKR-O. 0 
TOTAL-O.O 
ZPK-0.0 
2TIME-0.0 

RK I LL-R ‘DAMRAD/ 10.0 

R-R+1.0 

A-0.0 

00  151  J- UNARC 
THET  A=A*ALPHA 
ANGLE- THETA* 180. 0/P I 

AVSPO-F VSUB* ABSFl 1. 0-0. 9*THET A/ ( 2. 0*P I ) ) 

TO«TRAD*THETA/AVSPD 
A- A* 1.0 
8=0.0 

DO  150  K-1.M8T 
T1-B*DT 
T2-BT-T1 
T3-.T2-T0 
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CALL  GRNUNB (A  L ) 

DELTA*A1*SA 
CALL  GRNUMB ! A2 ) 

RTRUE=ABSF(50.0oA2»SR) 

ROBS-RrRU£+-A2*SR 
XOBS-“ROBS*SI NF! DELTA ) 

YOBS=  ROBS»COSF(UELTA)-RTRUE 
IF(r3)10, 15,20 

10  XSUB I ~  TRAD*  ! 1 .O-COSF !  AVSPD*T2/TRAD)L 
YSUBL*FVSUB»ri+TRAD*SI.NF!AVSPD»T2/TRAD) 

GQ  TO  25 

15  XSUBI^TRADMl.O-COSFITHETA)) 

YSUBI=FVSUB*T1+SINF( THETA) *TR  AO 
GO  TO  25 

20  XSUB I  *  TRAD* ( l.O-COSFi THETA) )  *FV SUB *T3*S I NF ( THET A ) 

YSUB  I.*FVSU6*T  1 *  TRADES  I NF  (  THETA  )*FV  SUB*  T3*C0SF!  THETA) 

25  DSUBI=SQRTF[ (XSUB I-XOBS) 4*2*1 YSUB! -YOBS)* *2) 

OPERS*DSUBI/S 

DRPERS*RKILL/S 

CALL  OCIP(OPERS,URPERS,PK> 

SUM=SUM*PK 

rOTAL=TOTAL+1.0 

B*B*1.0 

IFIPBLT- 1.0)26, 28, 26 
28  LFI63.0-C 135,35,30 

30  PRINT  31,  Tl, ANGLE, XSUB!., YSUBI.bSUUI.S.RKTlL, SPEED, PK 

31  FORMAT ( IH  .8F10.2,  F10.3) 

C-C+1.0 

GO  TO:  150 
35  PRINTi  52 

42  FORMAT  I 1H1.1H  ) 

PRINT  6 

C  =  S.O 
GO  TO  30 

26  |F(PBLT-2. 0)150, 27, 150 

27  I.F  I  71)150,37,  150 
37  ZPK*ZPK*PK 

7.T  I  MEALTIME*  1.0 
GO  TO  28 

150  CONTINUE 

151  CONTLNUE 

LF I PBIT-2.0 1162,160, 162 

160  4PK*Z;PK/Z  TIME 

I F I  6 3 • 0-C  1  159,  159,  161 
159  PRINT  42 
PRINT  6 
C*5.0 

161  PRINT  43. SPEED, RK I LL , S ,ZPK , ZT I  ME 

43  FORMAT!  1H  ,  6HSP£ED*F6 .2  ,.2X  ,  6HRK 1  LL  =  F  10 .2 , 2X ,  2HS*F  1 0. 2 , 2X.  6HZPKA  V- , 
1F6.3,.2X,6HT0TAL»F6.0) 

C*C*1.0 

162  AV.PK*SUM/TOTAL 
LFIPBIT-l. 0)171, 163, 163 

163  IF  1 63. 0-0170, 170, 171 

170  PRINT  42 
PRINT  6 
C».5.0 

171  PRINT  45, . SPEED, RK I LL,S,AVPK, TOTAL 

45  FORMAT!  1H  6HSPEED*F6.  2, 2X ,  6HRK I.LL*F;10. 2 , 2X, 2HS»F  10. 2, 2X, 5HAVPK-F6. 
13,2X,6HT0TAL»F8.0) 

C«C*l.O 

IFIPBIT- 1.0)152, 172, 172 

172  IFI63.0-C-2.0)  173,173,174 

173  PRINT  42 
PRINT  6 
C-5.0 
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GO  TO  152 

174  PRINT  6 
C=C+1.0 

152  CONTINUE 

IFI63.0-  0175,175,176 

175  PRINT  42 
05.0 

GO  TO  153 

176  PRINT  41 
C=C*1.0 

41  FORMAT  MHO,, -1H  ) 

153  CONTINUE 
NCASES=NCASES-l 

LF INCASES) 154, 154,2 

154  I  F  INSETS- 1)  50, 22, 50 
50  CALL  END  JOB 

END 
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l  1  0,  OUTPUT  OPTION  A 

COMPUTATION  OF  EFFECTIVE  KILL  PROBABILITY 


INPUT  PARAMETER  SET 

TURN  ANCLE  MULTLPLES  BUND  MULTIPLES  TURN  SU8  DAMAGE 

THETA  OF  THfeTA  TIME  OF  TIME  RADIUS  SPEED  RADIUS 

180.00  9.00  24.00  4.00  200.00  .40.00  10000.00 


AZIMUTH  RANGE 
ERROR  ERROR 

5.00  50.00 


BALLISTIC  TARGET. 

ERROR  ERROR 
50.00  50.00 


COMPUTATION  SETS 


SPEED* 

8.00 

RKILL* 

1000.00 

s* 

7Aa  7 1 

AVPK* 

1.000' 

/tOTAL- 

50. 

SPEED* 

8.00 

RKILL- 

2000.00 

s* 

70.71 

AVPK* 

1.000 

TOTAL- 

50. 

SPEED* 

8.00 

RKILL- 

3000.00 

s« 

70.71 

AVPK* 

1.000 

TOTAL* 

50. 

SPEED* 

8.00 

KKILL- 

4000.00 

S* 

70.71 

AVPK* 

1.000 

TOTAL- 

50. 

SPEED* 

8.00 

RKILL- 

5000.00 

s« 

7  0.7-1 

AVPK* 

i.ooo 

TOTAL* 

50. 

SPEED* 

8.00 

RKILL- 

6000. 00 

s- 

70.71 

AVPX- 

1.000 

TOTAL- 

50. 

SPEED- 

8.00 

RKLLL- 

7000.00 

s* 

70.71 

AVPK- 

1.000 

TOTAL* 

50. 

SPEED* 

8.00 

RKILL- 

6000.00 

s« 

70.71 

AVPK* 

1.000 

TOTAL- 

50. 

SPEED- 

8.00 

RKILL- 

9000. 00 

s- 

70.71 

AVPK® 

1.000 

TOTAL- 

50. 

SPEED* 

8.00 

RKLLL* 

10000.00 

s» 

70.71 

AVPK* 

1.000 

TOTAL- 

.50. 

SPEED- 

U.00 

RK  LU¬ 

1000.00, 

s- 

70.71 

AVPK* 

0.969" 

TOTAL- 

50. 

SPEED* 

16.  OO1 

RK  ILL'*' 

2000.00 

s- 

70.71 

AVPK* 

1 . 000 

TOTAL* 

50. 

SPEED- 

16.00' 

RKILL- 

3000.00 

s* 

70.71 

AVPK* 

1.000 

TOTAL- 

50. 

SPEED* 

16.00’ 

RKILL- 

4QOO.OO 

s- 

70.71 

AVPK* 

1.000 

TOTAL* 

50. 

SPEED- 

16.00 

KKLLL* 

5000. 00 

s- 

70.71 

AVPK* 

1.000 

TOTAL- 

50., 

SPEED* 

16.00 

RKLLL- 

6000. 00 

s* 

70.71 

AVPK* 

1.000 

TOTAL* 

50. 

SPEED* 

16.00 

RKILL* 

7000.00 

s* 

70.71 

AVPK* 

1.000 

TOTAL- 

'SO. 

SPEED- 

16.00 

RKILL* 

8000.00 

s* 

70.71 

AVPK* 

.1.000 

TOTAL* 

50. 

SPEED- 

16.00: 

RK1LL- 

9000. 00 

s- 

70.71. 

AVPK* 

1.000 

TOTAL- 

SO. 

SPCED* 

16.00 

RK1LL* 

10000.00 

S-, 

70.71 

AVPK* 

1.000 

TOTAL- 

50. 

SPEED- 

24.00 

KKILL- 

1000.00 

s- 

70.71 

AVPK* 

0.768 

TOTAL* 

50. 

SPEEO- 

24.00 

RKLLL- 

2000.00 

s- 

70.71 

AVPK* 

1.000 

TOTAL- 

50. 

SPEED- 

24.00 

RKILL* 

3000.00 

s« 

70.71 

AVPK« 

1.000 

TOTAL* 

50. 

SPEED- 

24.00 

RKILL- 

4000.00 

s« 

70.71 

AVPK* 

1.000 

TOTAL* 

50. 

SPEEO- 

24.00 

RKI.LL- 

5000.00 

s* 

70.71 

AVPK— 

1.000 

TOTAL- 

50., 

SPEEO- 

24.  GO 

RKILfc- 

6000.00 

s- 

70.71 

AVPK* 

1.000 

TOTAL* 

50. 

SPEED- 

24.00 

RK1LL- 

7000.00 

s- 

70.71 

AVPK*' 

1.000 

TOTAL* 

.50. 

SPF.ED- 

,24.00 

RKILL* 

8000.00 

s- 

70.71 

AVPK- 

1.000 

TOTAL* 

50. 

SPEED- 

24.00 

RKILL*' 

9000.00 

s- 

70.71 

AVPK* 

1.000 

TOTAL* 

50. 

SPEED* 

24.00 

RKILL- 

10000.00 

s- 

70.71 

AVPK- 

1.000 

TOTAL- 

50. 

SPEED- 

32.00 

RKILL* 

1000.00 

s- 

70.71 

AVPK* 

0.215 

TOTAL* 

50. 

SPEED* 

32.00 

RKLLL- 

2000.00 

s- 

70.71 

AVPK* 

1.000 

TOTAL* 

50. 

SPEEO- 

32.00 

RKILL- 

3000.00 

s- 

.70.71 

iVPK>. 

1=000 

TOTAL- 

5Q» 

SPEEO* 

32.00 

RKILL- 

4000.00 

s* 

70.71 

AVPK*- 

1.000 

TOTAL* 

50. 

SPEEO* 

32.00 

RKILL* 

5000.00 

s- 

70.71 

'AVPK- 

1.000 

TOTAL® 

50. 

SPEED- 

32.00 

RKILL - 

6000.00 

s- 

7C.71 

AVPK* 

1.000 

TOTAL- 

50. 

SPEED- 

32.00 

RKIU * 

7000.00 

s- 

70*71 

AVPK* 

J..0QO 

TOTAL- 

SO. 

SPEED- 

32.00 

RKILL* 

8000.00 

s- 

7-5  ,?i 

AVPK* 

,1.000 

1QTAL- 

50. 

SPEED- 

32.00 

KKILL- 

9000. 00 

s- 

AVPK- 

1.000 

TOTAL- 

50. 

SPEED* 

32.00 

RKILL- 

10000.00 

s® 

70.71 

AVPK- 

1.000 

TOTAL- 

50. 

SPEED- 

40.00 

KKLLt- 

1000. 00 

s* 

70.71 

AVPK- 

0.076 

TOTAL* 

50. 

SPEED* 

40.00 

HKILL- 

2000.00 

s- 

70.71 

AVPK* 

1.000 

TOTAL- 

50. 

spetd- 

40.00 

RKI.LL- 

3000. 00 

s- 

70.71 

AVPK- 

1.000 

TOTAL* 

50. 
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L  1  1.  OUTPUT  OPTION:  6 

COMfUTAIION  OP  EFFECTIVE  KILL  PRUB AB II. I  TV 


INPUT  PARAMETER'  SET 

fURN  ANGLE  MULTIPLES  BLIND  MUL-MPfeES  -TURN  SUB  DAMAGE  A/.  JMUTH  RANGE  BALLISTIC  TARGET 


THETA 

uo.oo 

OP  THETA 
9.00 

TIME  UP  TIKE  RADIUS 

24.00  4.00  200.00 

SPEED 

40.00 

RADIUS 

10000.00' 

ERROR 

5.00 

ERROR 

50.00 

ERROR 
50;  00 

ERROR 

50.00 

TIME 

THETA 

COMPUTATION  SETS. 

XSUB!  V  SUB  I  DSUQI  SIGMA 

RKILL 

SPEED 

PK 

0. 

G. 

-0. 

324.29 

6  39.  12 

70.71 

1000.00 

8.00 

1.000 

6.00 

0., 

-0. 

324.29 

318.01 

70.71 

1000.00 

8.00 

1.000 

12.00 

0. 

-0. 

324.29 

338.05 

70.71 

1000.00 

8.00 

1.000 

I  B .  00 

0. 

-0. 

324.29 

334.38 

70.71 

1000.00 

8.00 

1.000 

24.00 

0. 

-0. 

324.29 

360.87 

70.71 

1000.00 

8.00 

1. 000 

0. 

20.00 

71.6V 

302.78 

4  18.-13 

7.0.71 

1000.00 

8.00 

1.000 

6.00 

20.00 

43. -77 

307.67 

238.15 

70.71 

LOOb.CO 

•8.00 

1.000 

12.00 

20.00 

19.66 

314.49 

339.77 

70.71 

1000.00 

8.00 

1.000 

18.00 

20.00 

4.V4 

320.02 

732.39 

70.71 

1000.00 

8.00 

T.  .000 

24.00 

20.00 

-6. 

324.29 

302.97 

70.71 

1 000. 00 

8.00 

1.000 

0. 

40.00 

69.60 

280.48 

294.46 

70.71 

1000.00 

8.00 

i-.ooo 

6.00 

40.00 

39.49 

295.14 

404.01 

70.71 

1000. OU 

8.00 

1.000 

12.00 

40.00 

17.66 

306.64 

339.35 

70.71 

1000.00 

8.00 

I.  000 

18.00 

40.00 

4.43 

316.00 

347.29 

70.71 

1000.00 

,8.00 

1.000 

24.00 

40.00 

-0. 

324.29 

441.99 

70.71 

1000.00 

8.00 

1 . 000 

0. 

60 . 00 

62.2I.1 

266.05 

245.61 

70.71 

1000.00 

8.00 

1.000 

6.00 

60.00 

35. 2i 

283.74 

317.63 

70.71 

1000.00 

8.00 

1.000 

12.00 

60.00 

15.76 

298.76 

328.87 

70.71 

1000.00 

8.00 

1.000 

18.00 

60.00 

3.95 

311.98 

333.42 

•70*71 

1000.00 

8-00 

i.OQG 

24.00 

60.00 

-0. 

324.29 

343.40 

70.71 

1000.00 

8.00 

wooo 

0. 

80.00 

55.22 

251.42 

305.41 

70.71 

1000.00 

8.00 

1 . 000 

6.00 

80.00 

31.27 

272.25 

280.83 

70.71 

1000.00 

8.00 

1.000 

12.00 

30.00 

13.97 

290.85 

520.  1  1 

70.71 

1000.00 

8.00 

1.000 

18.00 

80.00 

3.50 

307.95 

339.72 

70.71 

1000.00 

8.00 

1.000 

24.00 

80.00 

-0. 

324.29 

34K32 

70. 71 

3000.00 

8.00 

1.000 

0. 

100.00 

48.62 

236. 6  L. 

292.89 

70.71 

1000. 00 

8.00 

1.000 

6.00 

100.00 

27.52 

260.69 

230.34 

70.71 

1000.00 

8.00 

,1.000 

12.00 

100.00 

12.28 

282.92 

526.53 

70.71 

1000.00 

8.00 

i.ooo 

18.00 

100.00 

3.08 

303.92 

449.74 

70.71 

loob.oo 

8.00 

l  .•  000 

24.00 

100.00 

~0., 

324.29 

340.16 

70.71 

1000.00 

8.00 

1.000 

0. 

120.00 

42.43 

221.62 

2,55 . 83 

70.71 

1000.00 

8.00 

1-000 

6.00 

120.00 

23.99 

249.05 

240.40 

70.71 

1000.00 

8.00 

1.000 

12.00 

120.00. 

10.-70 

274.97 

296.00 

70.71 

1000.00 

8.00 

l'.OOO 

18. UO 

120.00 

2.68 

299.88 

308.78 

70.71 

1000.00 

8.00 

1.000 

24.00 

120.00 

-0. 

324.29 

639.52 

70.71 

1000.00 

8.00 

1.000 

0., 

140.00 

36.65 

206.48 

376.16 

70.71 

1000.00 

3-00 

1.000 

6.00 

140.00 

20.71 

237.34 

352.89 

70;  71 

1000.00 

8.  CTO  — 

1.000 

12.00 

140.00 

9.2T 

267.00 

273.03 

70.71 

1000.00 

8.00 

1.000 

18.00 

140.00 

2.31 

295.85 

42,9.49 

70.71 

1000*00' 

8.00 

1.000 

24.00 

14C.00 

-0. 

324.29 

320.21 

70.71 

1000.00 

8.00 

t .  000 

0* 

160  *  00 

31.2,7 

191 . 18 

279.09 

70.  71 

i 000.00 

8.00 

1 . 000 

6.00 

160.00 

i7 . 66 

225.57 

274.06 

70.71 

looo.oo 

8.00 

i.ooo 

12.00 

160.00 

7.87 

259.00 

279.98 

70.71 

1000.00 

8.00 

1.000 

18.00 

160.00 

1.97 

291.81 

4  71.26 

70.71 

1000.00 

0.00 

l'.OOO 

24.00 

160.00 

-0. 

324.29 

326.22 

70.71 

1000.00 

8.00 

1.000 

0. 

180.00 

26.32 

175.74 

191.71 

70.71 

1000.00 

8.00 

l .  000 

6.00 

180.00 

14.85 

213.74 

251.57 

70.71 

1000.00 

8.00 

1.000 

12.00 

uo.oo 

6.62 

251.00 

314.99 

70.71 

1000.00 

8. bo 

1.000 

C-4 
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2. 

OUTPUT 

OPTION;  C 

COMPUTATION  I 

OF  EFFECTIVE  KILL  PROBABILITY 

TURN  ANGLE  MULTIPLE’S 

BL  I'NU 

MULTIPLES. 

THETA 

OE  THETA 

T  IME 

OF  TIME 

180.00 

9.00 

24.00  4.0< 

COMPUTATION  l 

TIME 

THETA 

x  sue  i 

YSUBI  I 

0. 

0. 

-0. 

324.29 

0. 

20.00 

71.69 

302.78 

0. 

40.00 

69.60 

280.48 

0. 

60.00 

62.21 

266.05 

0. 

80.00 

53.22 

251.42 

0. 

100.00 

48.62 

2  36.61 

0. 

120.00 

42.43 

221.62 

0. 

L40.00 

36.65 

206.48 

0. 

160.00 

31-27 

191.18 

0. 

180.00 

26.32 

173.74 

SPEEOm  8.00 

KKLLL* 

1000.00 

S*  70, 

SPEEO*  8.00 

RK  ILL* 

1000.00 

S*  70. 

TIME 

THETA 

XSU8L 

YSUBI 

o.. 

0. 

-0. 

324.29 

0. 

20.00 

71.69 

302 .76 

0. 

40.00 

69.60 

280. 48 

0. 

60.00 

62.21 

266.05 

0. 

80^00 

33.22 

251.42 

0. 

100.00 

46.62 

236.61 

0. 

120.00 

42.43 

221.62 

0. 

140.00 

36.65 

206.46 

0. 

160.00 

31.27 

191.18 

0. 

180.00 

26.32 

175.74 

SPEEO-  8.00 

RKILL* 

2000;.  00 

S*  70. 

SPEEO*  8.00 

RKILL* 

2000.00 

S-  70. 

TIME 

THETA 

XSUbl 

YSUBI 

0. 

0. 

-0. 

324.29 

0. 

20.00 

71.69 

302.78 

0. 

40.00 

69.60 

28C.4B 

0. 

60.00 

62.21 

266.05 

0. 

80.00 

55.22 

251.42 

0. 

100.00 

48.62 

236.61 

0. 

120.00 

42.43 

221.62 

0. 

140.00 

36.65 

206.48 

0. 

160.00 

31.27 

191.18 

0. 

180.00 

26.32 

173.74 

SPEED*  8.00 

RKILL* 

3000.00 

S*  70. 

SPEEO*  R.00 

RKILL- 

3000.00 

S-  70. 

TIME 

THETA 

xsubi 

YSUBI 

0.. 

0. 

-0. 

324.29 

0. 

20.00 

71.69 

302.78 

0. 

40.00 

69.60 

280.48 

0. 

60.00 

62.21 

266.05 

0. 

80.00 

55.22 

251.42 

0. 

100.00 

48.62 

236.61 

0. 

120.00 

42.43 

221.62 

0. 

140.00 

36.45 

206.48 

0. 

160.00 

31.27 

191.18 

INPUT  PARAMETER  SET 

TUAN  SUB  DAMAGE  AZIMUTH;  RANGE 

RADIUS  SPEED  RADIUS  ERROR  ERROR 

200.00  40.00  LOOOO.OO  5.00  50.00 


3U.26 

312.95 
285.76 
3SS.Q7 

270.55 
240.07 
214.85 

188.69 
536.80 
4  10.13 

1  ZPKAV1 
l  AVPK* 
SUB! 

358.55 
456.28 
343.04 

272.69 
259.64 
204.39 
2  34.25 
225.90 
163.68 
365.37 

1  ZPKAV 
l  AVPK- 
SU61 

345.26 

264.95 

337.55 
293.01 
216.81 
286. 15 
303.14 
194.60 
237.02 
156.02 

ZPKAV' 
AVPK* 
DSUBI 
366.99 

303.26 
342.65 
299.78 
366.92 

220.98 

201.98 
212.64 
297.47 


71 


SIGMA  RKIU 

70.71  1000, 

70.71  1000, 

70.71  1000, 

70.71  1000. 

70.71  1000, 

70.71  1000, 

70.71  1000. 

70.71  1000 

70.71  1000 

70-71  iooo. 

»  1.000  TOTAL* 
1.000  TOTAL* 
SIGMA  RKRL 

TO. 71  2000 

70.71  2000. 

70.71  2000, 

70.71  2000 

70.71  2000 

70.71  2000. 

70.71  2000. 

70.71  2000. 

70.71  2000. 

7U,?l  2000. 

•  1.000  TOTAL- 
1.000  TOTAL* 
SIGMA  RKILl 

70.71  3000. 

70.71  3000. 

70.71  3000. 

70.71  3000. 

70.71  3000. 

70.71  3000. 

70.71  3000. 

70.71  3000, 

70.71  3000, 

70.74  3000, 

»  1.000  TOTAL* 
1.000  TOTAL- 
SIGMA  RKILL 

70.71  4000 

70.71.  4000 

70.71  4000 

70.71  4000 

70.71  4000 

70.71  4000 

70.71  4000, 

70.71  4000 

70.71  4000 


SPEED 

,00  8.00 

,0U  8,00 

,00  8.00 

.00  8.0Q 

.CO  8.00 

,00  8.00 

00  6.00 

00  8.00 

00  8.00 

00  8.00 

10. 

50. 

SPEEO 

.00  8.00 

00  8.00 

,00  8.00 

,00  6.00 

.00  8.00 

.00  8.00 

,00  8.00 

,00  6.00 

.00  8.00 

.00  8.00 

10. 

50. 

SPEED 

.00  8.00 

.00  8.00 

.00  8.00 

.00  8.00 

.00  6.00 

.00  8.00 

00  8.00 

00  6.00 

00'  8.00 

.00  8.00 

10. 

50. 

SPEEO 

,00  6.00 

.00-  6.00 

00  6.00 

i.OO  8.00 

00  8.00 

,00  8.00 

,00  8.00 

,00  8.00 

,00  8.00 


PK 

1.000 
1.000 
1.000 
1.000 
1 .000 
1.000 
1.000 
1.000 
1 1 000 
1.000 


PK 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 


PK 

1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
1 . 000 
1.000 
1.000 
1.000 


PK 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 


BALLISTIC  TARGET 

ERROR  ERROR 
50.00  50.00 
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Y-axis 


Tactical  Situation: 


TB 

V 


R 

A 


blind  time 
sub's  speed 


N^  A  T  =  time  sub  begins  to  turn 

turn  radius 

Pj  A  0  =  turn  angle 


/  <xi' 


Yj)  represents  the  distri- 


Vbutions  of  points  in  a  lobe  of  the 
locus  of  an  evading  submarine  at 
the  end  of  the  blind  time. 


FIGURE  1 
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